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We consider a Potts model diluted by fully frustrated Ising spins. The model corresponds to 
a fully frustrated Potts model with variables having an integer absolute value and a sign. This 
model presents precursor phenomena of a glass transition in the high-temperature region. We show 
that the onset of these phenomena can be related to a thermodynamic transition. Furthermore this 
transition can be mapped onto a percolation transition. We numerically study the phase diagram 
in 2 dimensions (2D) for this model with frustration and without disorder and we compare it to the 
phase diagram of i) the model with frustration and disorder and of ii) the ferromagnetic model. 
Introducing a parameter that connects the three models, we generalize the exact expression of 
the ferromagnetic Potts transition temperature in 2D to the other cases. Finally, we estimate the 
dynamic critical exponents related to the Potts order parameter and to the energy. 
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I. INTRODUCTION 

The glass transition temperature T g for liquids is de- 
fined experimentally as the onset of calorimetric anoma- 
lies [jlj . It is usually understood that T g is not related to 
a thermodynamic transition Q , but to the slowing down 
of one or more degrees of freedom. This slowing down 
prevents the system to reach the equilibrium This 
question is of interest because almost all liquids can form 
glasses if cooled at high enough rate. Moreover, many 
other materials as for example polymers, microemulsions, 
granular material, vortex glasses, ionic conductors, col- 
loids, plastic glassy crystals and spin glasses (SGs) ||] 
show glassy properties. 

Even well above T g , where glassy systems actually can 
equilibrate, they show experimentally dynamic anomalies 
as precursor phenomena of the glass transition ||||-||. 
From a theoretical point of view, one of the open ques- 
tions is if these precursor phenomena are related to the 
thermodynamics of the system [ppP,^0[|, or if they are 
not, like the glass transition occurring at lower temper- 
ature T g . In particular, for the Ising SG || such a rela- 
tion between precursor phenomena and a thermodynamic 
free energy essential singularity has been shown P, pr| , ^2| . 
Indeed, in this model there is a dynamic anomaly at a 
temperature T* . Above T* the relaxation processes have 
an exponential behavior, while below T* they have a 
non-exponential behavior. Theoretical and numerical ev- 
idences show that T* coincides with the Griffiths temper- 
ature T c |Tl| . This T c is the transition temperature that 
the model would have if the frustration due to disorder is 
removed. A way to remove the frustration is, for example, 



to substitute every antiferromagnetic interaction with a 
ferromagnetic interaction. In general, removing the frus- 
tration, the model will have ferromagnetic regions and 
antiferromagnetic regions and T c will be the transition 
temperature of the unfrustrated model. To be more pre- 
cise, the free energy of the Ising SG in external field has 
a singularity which disappears in the limit of zero exter- 
nal field and which occurs at a temperature that goes to 
T c in the same limit. This singularity is present only in 
disordered systems. The relation T* = T c has also been 
shown to be valid numerically in 2 dimensions (2D) for a 
model with Potts variables that generalizes the Ising SG 
fUjfl . This generalized model will be considered in the 
following and we will refer to it as the Potts SG. 

Until now we talked about glass dynamics in systems 
with disorder, but numerical simulations show that it is 
possible to observe glassy behavior with precursor phe- 
nomena for spin systems without disorder but with frus- 
tration piyiq] . These systems, due to the lack of disor- 
der, are more suitable for a theoretical approach p6| , p7[ . 
In particular, one can try to answer the question about 
the relation between the precursor phenomena and the 
thermodynamics of the system. 

In Ref. JTJ] it was considered a simple case: The fully 
frustrated (FF) Ising model |L9| where ferromagnetic and 
antiferromagnetic interactions are ordered in such a way 
that any lattice cell has an odd number of antiferromag- 
netic interactions (i.e. is frustrated). It is shown (by 
simulations in 2D and 3D) that the onset T* of non- 
exponential relaxation processes is related to a random- 
bond percolation transition |Q. As for the Ising SG, it is 
possible to generalize the FF Ising model to a FF model 
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with Potts variables Q: The Potts FF model that will 
be described in details in the following. In Ref. Jl3| the 
dynamics of the Potts FF model was compared with the 
dynamics of the Potts SG model and some anticipations 
on the relation between precursor phenomena and ther- 
modynamics were given. In this work we study in details 
the thermodynamics of the Potts FF model, showing that 
T* corresponds to the thermodynamic transition temper- 
ature T p of the Potts variables. It is important to note 
that in any FF model the Griffiths temperature T c cannot 
be defined for the lack of disorder. Therefore T c cannot 
play any role in these cases. 

The comparison of the results presented here with the 
analogous study of the Potts SG model with disorder 
fill allows to give insight on the role of disorder. To this 
goal we introduce a formal parameter < X < 1 that 
connects both models. For X — we have the ferromag- 
netic Potts model |2(| (without frustration and without 
disorder). For < X < 1 we have the disordered and 
frustrated Potts model. In particular, for X — 0.5 we 
have the Potts SG model. For X = 1 we have the Potts 
FF model (without disorder). 

From Ref. [jl8) and from the present work, it is pos- 
sible to show that for X = 0.5 and X = 1 there are 
two thermodynamics transitions. The lower transition is 
an Ising SG or a FF Ising transition, respectively. The 
upper transition at T p is in the universality class of a fer- 
romagnetic Potts transition. Furthermore we show that 
T p corresponds to a percolation temperature. Moreover, 
we show how it is possible to generalize the exact expres- 
sion of T p for the model with X = in 2D |20) to the 
cases X = 0.5 and X = 1. 

The organization of the paper is the following. In 
Sec. II we introduce the model and the known results for 
X = (ferromagnetic case), for X = 0.5 (disordered and 
frustrated case) and for X = 1 (ordered and frustrated 
case). In Sec. Ill we introduce the cluster formalism used 
to map the upper thermodynamic transition at T p onto 
a percolation transition. In Sec. IV we present the phase 
diagram in 2D for X = 1 as result of Monte Carlo (MC) 
simulations and we compare it with the cases X = and 
X = 0.5. In Sec. V we use the spin-flip MC dynamics to 
study the dynamic critical exponent and the temperature 
T*, onset of stretched exponentials. In Sec. VI we give 
the summary and conclusions. 



II. THE MODEL 

Structural glasses, such as dense molecular glasses, 
plastic crystal, or orto-therphenyl at low temperature, 
can be modeled to first approximation as systems with 
orientational degrees of freedom frustrated by geomet- 
rical hindrance between non-spherical molecules. For 
this reason we will consider the lattice model intro- 
duced in Ref. j2l), where the orientational degrees of 
freedom are represented by Potts variables pQ] with s 



states (<Ji = l,...,s) and the frustration is modeled 
by means of ferro/antiferromagnetically interacting Ising 
spins (Si — ±1), coupled to the Potts variables. 
The model is defined by the Hamiltonian 

II .{S,.n,.i , _ | = - S J /^^^(ejjSjSj + 1) (1) 

where the sum is extended over all the nearest neigh- 
bor sites, J is the strength of interaction, eij = ±1 
is a quenched variable that represents the sign of the 
ferro/antiferromagnetic interaction and 8 n , m = 0, 1 is a 
Kronecker delta. To emphasize that the Ising and the 
Potts variables arc interdependent, we can rewrite the 
Hamiltonian in Eq.([l]) as 

H s {n,e id } = 12,./)^, V -. r (2) 

where Ti = SiUi — ±1, ±2, . . . , ±s is a variables with 2s 
states and the frustration now is explicitly on the new 
variable. Any Ti has an absolute value Ui and a sign Si. 

The model depends on the interaction configuration 
{eij}. Possible choices are the following, (i) If all e,,j = 1 
(i.e. all interactions are ferromagnetic), the Eq.fEr) is the 
Hamiltonian of the ferromagnetic Potts model EG] with 
variables Ti with an even (2s) number of states. It shows a 
thermodynamic transition at T c (s) whose order depends 
on s. (ii) If £i j are quenched random variables, the model 
corresponds to the Potts SG. This model is a generaliza- 
tion of the Ising SG model that is recovered for s = 1. It 
shows two thermodynamic transitions |1S| ]. The lower is 
a SG transition at Tsg(s)- The upper is a Potts transi- 
tion at Tp(s) > Tsg(s). The transition at T p (s) is in the 
universality class of a ferromagnetic s-states Potts model. 
Another relevant temperature for this model is the tem- 
perature T c (s) defined for the previous case (the ferro- 
magnetic 2s-states Potts model). Indeed, it is possible 
to show that for finite external field a free energy (Grif- 
fiths) singularity arises . In the limit of external field 
going to zero, the temperature at which this singularity 
occurs goes to T c (s) and the singularity vanishes. T c (s) 
is the Griffiths temperature for this model. Furthermore, 
numerical simulations |l3f| show that T c (s) > T p (s) corre- 
sponds to the onset T*(s) of non-exponential correlation 
functions for the Ising spins Si. This result generalizes 
what happens in the Ising SG (s = 1 case), where at 
the Griffiths temperature T c (l) non-exponential corre- 
lation function are seen jig] . It is worth to note that 
the Ising spins are critical at Tsg(s), that is well below 
T c (s) = T*(s). Moreover, note that the relevant Grif- 
fiths temperature for this model is T c (s) of the variables 
Ti = SiCTi and not the Griffiths temperature of the vari- 
ables Si alone, that in our notation is T c (l). (Hi) If there 
is an odd number of e,-,j = — 1 for each elementary cell (as 
in Fig. |l|), the system is fully frustrated (FF). It means 
that at least one interaction per cell is not satisfied, i.e. 
the relative energy contribution is instead of —2s J (as 
for the edges b, d and / in Fig. |l|). The model is called 
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Potts FF and is a generalization of the FF Ising models 
|l9f that is recovered for s = 1 . It has frustration but no 
disorder. For any integer s > 1 the model has in 2D a 
second-order phase transition at T = 0. Non-exponential 
correlation functions are reported below a finite tempera- 
ture T*(s) for s = 2, 1 and 1/2 (the latter case is defined 
in Sec. Ill) (14||k|. This dynamic anomaly cannot be re- 
lated, as in the previous case, to the Griffiths tempera- 
ture. Indeed, the Griffiths temperature is not defined in 
FF models for the lack of disorder. The aim of this work 
is to study the phase diagram of the FF model as func- 
tion of s and to show that T*(s) corresponds to a Potts 
transition. Furthermore, using a percolation approach, it 
is possible to show that the Potts transition corresponds 
to a percolation transition, defined also for non-integer 
values of s. 

^ ., f b >5 c Vl \ 

} * \ \ } 

FIG. 1. Example of Potts fully frustrated model 
on a square lattice: on each vertex there is variable 
Ti = Sicn — ±1, ±2 . . . ± s with s = 4 in the figure. Here 
we represent the sign (Si) of each t; by an open or a full 
dot (respectively positive and negative, for example) and its 
orientational state (oi) by an arrow pointing in 4 different di- 
rections. Ferromagnetic (antiferromagnetic) interactions are 
represented by full (dotted) lines. 

The previous three cases can be generalized in a class 
of Hamiltonians. Consider, as model parameter, the den- 
sity X of frustrated elementary cells. For X — there 
are no frustrated cells (like in a ferromagnet), while for 
X = 1 every cell is frustrated (FF case) . For < X < 1 it 
is possible to partition the lattice in two non-overlapping 
sub-sets U(X) and F(X) of unfrustrated cells and frus- 
trated cells, respectively. Therefore the Eq.(|J) can be 
written as 

H s , x { n } = H s { n , e id (X)} = -2sJ £ 6 Ti , T . 

(i,j)eu(x) 

-2s J Yl S ^n,r } • (3) 

(i,j)£F(X) 

In this way the ferromagnetic Potts model with an even 
number of states (case i) is recovered for X — 0, since 



F(X — 0) is empty; The Potts SG model (case ii) for 
X = 0.5; The Potts FF model (case Hi) for X = 1, since 
U(X = 1) is empty; Intermediate disordered models are 
obtained for other values of < X < 1. 

To make the comparison between the X = 0.5 and 
X = 1 model, we will follow the line of Ref. jl8| where 
the phase diagram of the former has been studied. In 
particular, for X = 1 we consider a square lattice with 
one 6ij = — 1 per cell as shown in Fig.[l]. 

III. PERCOLATION 

We now introduce the percolation map following the 
Ref. It is a generalization to X > of Fortuin- 

Kasteleyn p2[ clusters formalism, defined originally only 
for X = 0. In Ref. |^l[ it is shown that the partition func- 
tion Z s ^x of the Hamiltonian in Eq.(^) can be written in 
terms of bond configurations C 



Z s . x = ]T e-^-^'A^ = ^ W S , X {C) 

{Ti} C 



-H S:X {r t }/(k B T) = w 



(4) 



where ks is the Boltzmann constant, W s ,x{C) = if C 
includes any frustrated loop (defined below), otherwise 



W S , X (C) =pW(l-p)^(2s) 



\A\, 



)N(C) 



(5) 



where p = 1 — exp[— 2s J/(fcsT)] is the probability to 
place a bond between two nearest neighbor sites, N(C) 
is the number of clusters (defined as maximal sets of con- 
nected bonds) in the configuration C, |C| is the number 
of bonds and |C| + 1 A\ is the total number of interactions. 
A loop of bonds is called frustrated if the product of all 
the signs tij of the interactions along it is equal to — 1. 
For a frustrated loop there is no {r^} configuration able 
to minimize the energy of all the interactions along it. An 
example of such a frustrated loop is the one composed by 
the edges a, e, g, f in Figjl] or a loop composed by the 
external edges of any three adjacent elementary cells in 
a FF lattice (or of any odd number of adjacent cells). 

Note that for X = there are no frustrated loops and 
the original Fortuin-Kasteleyn cluster definition is recov- 
ered. Since the above cluster definition holds for any s 
and X, for any model described by the general Hamilto- 
nian in Eq.(||) it is possible to define a percolation tem- 
perature for these clusters. 

Moreover, the right-hand side expression of Eq.([|) is 
well defined even for non-integer values of s. In these 
cases it does not define a Hamiltonian model, but it still 
defines a percolation model. In particular, the s = 1/2 
case for X > is the frustrated percolation model |TJ] . 

For X — 0, in the general case in which we consider 
Potts variables Ti with a generic number q of states in- 
stead of an even (2s) number of states in Eqs.(|3||5|), it 
is possible to show |2(| that the percolation tempera- 
ture T p (q, X — 0) of the above defined clusters coincides 
with the ferromagnetic g-state Potts critical temperature 
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T c (q). In particular, in 2D it is possible to prove (2(j] the 
relation 



k B T c (q) 



1 



qj ln(l + ^/g) 



(6) 



For X = 0.5 the same kind of relation between the 
percolation temperature T p (q,X = 0.5) and the Potts 
transition temperature for q = 2s in 2D has been ex- 
tended ]lS| . Furthermore, it was shown numerically that 
it is possible to generalize the Eq.(^) using a fitting pa- 
rameter a{X) Jig]. The resulting relation is 



k B T p (q,X) 
a{X)qJ 



1 



In 



1 + V«W9 



(7) 



with a(X = 0.5) = 0.800 ± 0.003 Q. Note that 
T p (q,X = 0.5) represents a percolation temperature for 
an y 9 > and also a Potts transition temperature for 
even integer values q = 2s. Furthermore, for < X < 1 
and q = 2s the Eq.(f3|) for T c (2s) gives by definition the 
Griffiths temperature of the disordered model. 



IV. NUMERICAL RESULTS IN 2D FOR X = 1 

To study the phase diagram of the model with X = 1 
(Potts FF) we have simulated it for s = 2, 7, 20, 50 
on square FF lattices with periodic boundary conditions 
and linear size L from 10 to 80. At low temperature all 
the Potts variables tend to order ferromagnetically wher- 
ever the Ising spins satisfy the ferro/antiferromagnetic 
interactions (i.e. S eij s it Sj — !)• I n particular, at T = 
the system in 2D has a second-order phase transition of 
the Ising spins, like in the FF Ising model (T^|. We will 
show that the interplay between Potts variables and Ising 
spins affects the phase diagram at finite temperature. As 
in Ref . |TJ] , to study the finite temperature range we can 
use an efficient cluster dynamics with an annealing pro- 
cedure. We define a MC step as an update of the whole 
system. At each temperature we average the data over 
10 4 MC steps, discarding the first 5 x 10 3 MC steps. 

In our systematic analysis we calculated for each s the 
Binder's parameter p4] for the energy density E defined 
as 



V = 1 



(E 4 



3(E 2 )< 



(8) 



(angular brackets denote the thermal average). This 
quantity allows to distinguish between first-order and 
second-order phase transitions. Indeed, in a second-order 
phase transition for L — > oo it is = 2/3 for all tempera- 
tures, while in a first-order phase transition V has a well 
pronounced minimum near the transition temperature. 
The thermodynamic of the model is studied by means of 
the Potts order parameter 



M 



s maxi(Mi) 



(9) 



(where i = 1, . . . s, M t is the density of Potts spins in the 
i-th state), the susceptibility 



X 



l (M 2 ) - (My 



N 



(10) 



(where N is the total number of Potts spins) and the 
specific heat 



C H = 



1 (E 2 )-(Ef 
k B T 2 N 



(11) 



To study the Fortuin-Kasteleyn percolation we calcu- 
lated the percolation probability per spin P = I — mi and 
the mean cluster size S = rn-i (where m n = ^ fc k n nk is 
the n-th moment of the distribution of density of clus- 
ters with size k). 



A. The second-order transition for s = 2 

For X — 1 and s — 2 the Binder's parameter V goes 
to the constant value 2/3 for all the temperatures as L 
increases (see Fig. |^) revealing a second-order phase tran- 
sition. The transition temperature T s in the thermody- 
namic limit can be estimated, together with all the crit- 
ical exponents, using the standard scaling analysis pq| 
for second-order phase transition, where by definition of 
critical exponents (3, 7 and a is £ ~ \T — T s \~ u be- 
ing £ the correlation length, M - \T - T s f - C"^, 
X ~ \T - T s |-T ~ C„ ~ |T - T s \- a ~ i a l\ for 

which we expect 



M~L-Pi»f M ((T-T s )L 1 ' 1 ') , 



(12) 



(and analogous scaling laws for \ and Cr) where /m(x) 
is an universal function of the dimensionless variable x. 
The values at which the scaling laws are verified give the 
estimates of critical exponents and of T s . 

In Figs. H we show the large-sizes data collapses us- 
ing the set of Ising critical exponents [^6| and leaving 
only T s as free parameter, giving an estimate ksTs/ J — 
2.73 ±0.03. 

The study of the percolation quantities shows a smooth 
behavior of P and a cusp in S increasing with L. There- 
fore one can make the ansatz that the percolation transi- 
tion is of second-order. Let us define a Fortuin-Kasteleyn 
percolation temperature T p and a set of percolation crit- 
ical exponents v p , (3 P and j p by means of the relations 
£ p ~ \T — T p \~ l,p where £ p is the connectedness length 
of the clusters (i.e. the typical linear cluster size) and 

P ~ \T - T p \Pp - i P Pvlv \ S ~ |T - T p |-> - £ p lvlVv . 
Applying the standard scaling analysis for percolation 
p7[ , one obtains in a consistent way the results sum- 
marized in Fig. ^ using the corresponding thermody- 
namic critical exponents for the 2D Ising model. The 
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numerical estimate for the percolation temperature is 
kBTp/J = 2.73 ± 0.03 coincident with the estimates of 
T s . Therefore the percolation transition and the Potts 
transition occur at the same temperature and with the 
same set of critical exponents, i.e. they coincide. 



k,Tp/J = 2.73±0.03 



> 
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0.66 




k B T/J 

FIG. 2. X = 1 and s — 2: Binder's parameter V vs. T for 
the lattice sizes L listed in the figure. Inset: Enlarged view. 
Where not shown, errors are smaller than symbols size. Lines 
are only guides for the eyes. For increasing L, V goes to 2/3 
for all temperatures showing a second-order phase transition. 
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FIG. 3. X = 1 and s = 2: Collapse of M, Ch and x 
data in a), b) and c) panel, respectively, for the Ising criti- 
cal exponents (y, f3, a, 7). Each collapse gives an indepen- 
dent estimate of the critical temperature T s (indicated in each 
panel) . 
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panel a) and b), respectively, for the Ising critical exponents. 
Each collapse gives an estimate of the percolation tempera- 
ture k B T p /J = 2.73 ± 0.03. 



B. The first-order transition for s = 7, 20 and 50 

For X = 1 and s = 7, 20 and 50 we have considered 
systems with L from 10 to 50 lattice steps with periodic 
boundary conditions. On the base of the mean field re- 
sults [{l7| for the Potts FF model and the knowledge of 
the Potts model |£^| , we expect that the order of transi- 
tion changes for s > 4. In fact, the thermodynamic or- 
der parameter M becomes more and more discontinuous 
as s increases. At the same time the percolation order 
parameter P develops a better and better pronounced 
discontinuity. 

In particular, the study of the Binder's parameter V 
reveals that the model for s = 7, 20 and 50 has a first- 
order phase transition, since, for each considered lattice 
size L between 10 and 50, V has a non vanishing min- 
imum, as shown in Fig. [| for s = 7. In this cases the 
estimates of infinite size transition temperatures T s (s) 
and T p (s), for the thermodynamic and the percolation 
transition, respectively, can be done through the relation 



Tmax(L) - Tmax(oo) — U 



(13) 



where D = 2 is the Euclidean dimension, T m ax(i) is 
the finite-size temperature of the maximum of Ch and 
S, respectively, and T m ax(oo) is the transition temper- 
ature in the thermodynamic limit. Therefore T s (s) and 
T p (s) can be evaluated by linear fits on a Log-Log scale 
with one free parameter. The results are summarized in 
Tab. I. For any s, T s (s) and T p (s) are consistent within 
the numerical error. 
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FIG. 5. X = 1 and s = 7: Binder's parameter 1/ as in 
Fig. 2. V has a well denned minimum for any size, showing a 
first-order phase transition. 



C. Phase diagram 

The numerical results in Sees. IV. A and IV. B give rise 
to the phase diagram for X = 1. It is qualitatively sim- 
ilar to the phase diagram for X = 0.5 In Fig.^ 
we compare both of them with the phase diagrams for 
X = 0. 

For X — (Potts model with 2s-states variables) 
for any integer s there is a paramagnetic/ferromagnetic 
phase transition at finite temperature T c (s). The tran- 
sition is of second-order for 2s < 4 and of first-order 
for 2s > 4 [§(|. T c (s) is given by Eq.(|) with q = 2s. 
It coincides with the Fortuin-Kasteleyn percolation tem- 
perature. The percolation temperature is defined even 
for non-integer values of s. 

For X — 0.5 (Potts SG model with frustration induced 
by disorder and with 2s-states variables) there are two 
phase transitions for integer values of s. The lower tran- 
sition is a SG transition. It is supposed at T$g = 
for any s in 2D and at Tsg(s) > in higher dimen- 
sions as in the Ising SG pq |. The high-temperature 
transition is a s-states Potts ferromagnetic transition oc- 
curring at T p (s,X) given by Eq.(^) with q = 2s and 
a(X = 0.5) = 0.800 ± 0.003 jjjf. It is a second-order 
transition for s < 4 and a first-order transition for s > 4. 
T p (s,X) corresponds also to the Fortuin-Kasteleyn per- 
colation temperature, that is defined even for non-integer 
s |l8| . This model is disordered and its Griffiths temper- 
ature is by definition the transition temperature of the 
2s-states Potts model T c (s), given by the Eq.(||) with 
q = 2s. It is possible to see that T c (s) is numerically 
consistent with the onset of non-exponential relaxations 
at T = T*(s) for s = 2 in 2D jlf] and for s = 1 (Ising 



SG) in 2D and 3D Jl|g§. 

For X — I (Potts FF model with frustration and no 
disorder and with 2s-states variables), considered here, 
there are two transitions for integer values of s, as well 
as for X = 0.5. The lower transition is at Tpp = in 
2D and at Tff(s) > in higher dimensions, as for the 
FF Ising model (s = 1) [[19 32|. As seen in this section, 



in this case the upper transition is a s-states Potts ferro- 
magnetic transition at T p (s, X = 1). As shown in Fig. [|, 
Eq.(0) describes well T p (s,X = I), using a fit param- 
eter a(X = I) = 0.690 ± 0.003. Furthermore we have 
shown that it coincides with a Fortuin-Kasteleyn perco- 
lation transition. As consequence the Fortuin-Kasteleyn 
clusters represent the regions of correlated Potts vari- 
ables. This means that the cluster's characteristic linear 
size is equal to the correlation length of the Potts vari- 
ables. Analysis for s = 1/2 (frustrated percolation) and 
s = I (FF Ising model) in 2D and 3D are given in Ref. 
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FIG. 6. Numerical phase diagram in 2D for X = (Potts 
model: T c ), X = 0.5 (Potts SG model: T p PSG ) and X = 1 
(Potts FF model: T^ FF ). The data are fitted with Eq.(§) 
with fit parameter a(X) shown in the inset. Errors are smaller 
than symbols size. 

The main difference between X — 1 and X — 0.5 cases 
is that the Griffiths temperature is defined in the latter 
(disordered model), but not in the former (for the lack of 
disorder). The important consequence of this fact is that 
for X = 1 the dynamic anomalies are present only below 
T p (s, X = 1) gfj, while for X = 0.5 they are present also 
above T p (s,X = 0.5) and below T c (s) ||. 

Note that all the Potts-percolation transition temper- 
atures for X = 0, 0.5, 1 can be described by the same 
form in Eq.(^|) as function of s with a A-dependent pa- 
rameter a(X). This parameter has a regular behavior as 
function of X (see the inset in FigJ|). 
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FIG. 7. X = 1 and s — 2: The logarithm of correlation 
time tq (see text) at the finite-size transition temperature 
T P (L) p3 as function of the logarithm of size L, for the cor- 
relation function of the Potts order parameter M (circles) and 
of the energy density E (squares), for L = 20, 24, 30, 40, 50. 
The slopes give the exponents zm and ze- 



V. DYNAMIC CRITICAL EXPONENT 



As shown in Sec. IV.A, at T p (s = 2, X = 1) the Potts 
variables have a second-order phase transition. There- 
fore their correlation length diverges and their dynamics 
slows down. A measure of the slowing down for any ob- 
servables A is given by the dynamic critical exponent z 
defined by ta(T p (L),L) ~ L z . Here, omitting for sake 
of simplicity the dependence on s and A, T p {L) is the 
transition temperature for the system with finite size L 
p3| and ta [T, L) is the correlation time at temperature 
T and size L associated to the correlation function for A 

f (t T\ — (Mt,T)A(0,T)) - (A(T))i 

™ ' ' (A(0, T) 2 ) — (A(T)) 2 ' { 1 

where t is the time. More then one definition of ta is 
possible for any A, but all of them, even if numerically 
different, have the same qualitatively behavior (l2). In 
particular, facing the difficulty that the greater size the 
greater ta{T p (L),L), we define tm and te for Potts or- 
der parameter M and for energy density E as the time 
t (in unit of MC steps) at which /m(to> T p (L)) = 0.4 
and /e(tq, T p {L)) = 0.3, respectively. The data for sizes 
L = 20, 24, 30, 40, 50 are shown in Figj^ and the dynamic 
critical exponents are estimated by linear fits on Log-Log 
plot as zm — 0.9 ± 0.1 and ze — 1.1 ± 0.1, respectively 





t(MC steps) + 1 

FIG. 8. X — 1 and s — 2: Correlation function of the 
Potts order parameter M in the thermodynamic limit (see 
text). For clarity we show only some of the recorded data 
for the simulated temperatures. Symbols are results of sim- 
ulations, solid lines are fits with the form in Eq.(|l5|), dash 
lines (on this scale indistinguishable from solid lines) with 
stretched exponential form and dotted lines with exponen- 
tial form. Where not shown, the errors are smaller that the 
symbols size. Temperature is measured in J/fcs- 

To study the behavior of /m and Je in the thermody- 
namic limit we have extrapolated the data in the infinite 
size limit following the procedure suggested in Ref . |36[ . 
It consists in plotting at any t the generic /a(£, T, L) for 
finite size L versus 1/L and in extrapolating for 1/L — ► 0. 
We consider the temperature range 2.65 < ksT/ J < 3.25 
pEfl . The results are shown in Figs. || and ^|. To check 
the form of /m and /e, we have fitted the data with 
three different plausible functions: a) Simple exponen- 
tial /o-exp(— t/r); b) Kolraush- Williams- Watts stretched 
exponential fo ■ exp[— (t/r) 13 ]; c) Ogielski form 

e -(V^)' 3 

f(t,T) = f —^— (15) 

that is a combination of a form b) with a power law. In 
the previous functions, fo, r, x and j3 are T-dependent 
parameters. Note that for the form c) it is possible to es- 
timate fo and x separately from (3 and r, since the first 
two describe the short-time behavior while the second 
two the long-time regime. 

The Ogielski form turns out to describe very well fu 
and /b as shown in Figs. || and [| The fitting parameters 
are presented in Fig.^J and [l0|. The forms b) and c) give 
always compatible estimates of /3 exponent, consistently 
with the very low values of x (that is approximately the 
slope of the function at t = in Figs.|| and ||). 

From Fig. [lO] for the parameter (3 of /m, we see that 
it is j3 = 1 above the Potts transition at ksT p /J = 
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2.73 ± 0.03 and (3 < 1 below T p . Therefore there is a 
dynamic transition between a high-temperature exponen- 
tial behavior to a low-temperature stretched exponential 
behavior of /m- This result shows the presence of a com- 
plex dynamics below T p , consistently with the analysis of 
non-linear susceptibility correlation function of the Ising 
spins It is important to note that, while the Potts 
variables have a transition at T p , the Ising spins have no 



transition at T p and, in principle, 
lies are expected for them. 



no dynamical anoma- 




1 10 10 10 

t(MC steps) + 1 

FIG. 9. X — 1 and s = 2: Correlation function of the en- 
ergy density E in the thermodynamic limit as in Fig. 8. Inset: 
The correlation times estimated by exponential fit (circles) 
and by Eq. ( |l6| (squares). The arrow shows the numerical 
estimate of T p . 

On the other hand, the data for fg show that the long- 
time behavior is well described by an exponential func- 
tion. In particular, the form c) with /3 = 1 and a non-zero 
x fits very well also the data for the lowest temperature 
considered here. Note that, in principle, one can expect a 
dynamical anomaly even for the energy density E, since 
it depends explicitly on the critical Potts variables. 

In the fitting forms used, the correlation time r is de- 
fined as a fitting parameter. Another possible definition 
is the following (integral correlation time): 



Tint, A 



t r ) = JSo 5 + E/A(*,r) 



(16) 



Due to the divergence of the correlation time, the defi- 
nition in Eq.(|l6|) does not converge near the transition 
temperature. Therefore it is not shown in Figs. ^ and [l^ 
for the temperature with the largest r. However, where 
it converges, all the different estimates of the correla- 
tion time are numerically consistent as show in Fig.^ 
and uu. Note that, even if the data are extrapolated 



to the thermodynamic limit, finite size effects are still 
present. Indeed, the critical temperature estimated from 
r measurements, within our temperature mesh, is be- 
tween 2.75 J I ks and 2.85J/fcs, that is at variance with 
the estimate of T p . 
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FIG. 10. X = 1 and s = 2: Fit parameters for /m data 
with form in Eq. ( |l5| ) (circles) and for stretched exponential 
form (triangles). In the lower panel we show also the integral 
correlation times T; nt estimated with the Eq.(|l6|) (squares). 
The arrow and the vertical lines show T p and the numerical 
indetermination on it, respectively. 



VI. SUMMARY AND CONCLUSIONS 

We have studied a schematic model for glasses where 
frustrated orientational degrees of freedom, associated 
to 2s-states Potts variables Tj, induce complex dynam- 
ics. In real systems the frustration can be associated, for 
example, to steric hindrance of non spherical molecules 
in structural glasses. Each t% has an absolute value 
(o"j = 1, . . . , s) and a sign {Si = ±1). The frustration is 
over the signs and can be with disorder or without disor- 
der, depending on a model parameter X. For < X < 1 
there is disorder, while for X = 1 the model is (fully) 
frustrated without disorder. For X = there is no frus- 
tration, and the model recovers the 2s-states ferromag- 
netic Potts model. 

For < X < 1 the model has two thermodynamic 
transitions. The high-temperature transition at T p (s,X) 
for the model with 2s states is in the universality class 
of the s-states ferromagnetic Potts model. Therefore the 
fluctuations (x and Cjj) of the orientational degrees of 
freedom diverge at T p (s,X) (for s < 4), as well as the 
correlation times of quantities depending on them (like 
the Potts order parameter M or the energy density E). 
For them we estimated the dynamic critical exponents 
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Zm and ze- The low-temperature transition is a spin 
glass transition for X = 0.5, or a fully frustrated transi- 
tion for X = 1, and marks the ordering transition of the 
signs Si (Ising variables). 

The diverging fluctuations at upper and lower transi- 
tion temperatures are expected to be experimentally ob- 
servable only using specific probes that couple with them. 
Examples of such probes could be those associated to di- 
electric measurements in supercooled-liquids and plastic 
glassy crystals or to electron spin resonance spectroscopy 
measurements |^^]. 

For X = 1 and s = 2 the model shows a complex dy- 
namics in correspondence of T p . In particular, the cor- 
relation function for M and the correlation function for 
the signs Si change their behaviors. They have an expo- 
nential behavior above T p and a non-exponential behav- 
ior below T p . This behavior is expected at least for any 
s < 4, because it is associated to the free energy singular- 
ity occurring at the second-order transition of the Potts 
variables. Therefore the complex dynamics corresponds 
to a real thermodynamic transition. 

For X = 0.5 and s = 2 the onset of this dynamic 
anomaly is shifted to a higher temperature, above T p . It 
occurs in correspondence of the Griffiths essential singu- 
larity of the free energy. This singularity is not defined 
in the X = 1 case. For zero external field it goes to 
the transition temperature of the ferromagnetic 2s-states 
Potts model T c (2s) and it vanishes. Therefore the onset 
of complex behavior does not correspond to a real ther- 
modynamic transition. 

We have shown that the Potts transition at T p (s,X) 
for any s and X considered here, coincides with a per- 
colation transition. It is not worthless to note that the 
dynamic transition at T p (s, X = 1) persists also for s = 1 
and s = 1/2 Jl4|]. In these cases T p (s,X = 1) does not 
correspond to a thermodynamic transition, but only to a 
percolation transition. Therefore for s < 2 the dynamic 
anomaly is no more related to a thermodynamic transi- 
tion, but to a percolation transition in the real space |l4j . 
This result could be in some relation with experimental 
results on microemulsions 

Finally, we have shown that it is possible to general- 
ize the exact relation for the ferromagnetic Potts transi- 
tion temperature T c (s) in 2D to a transition temperature 
T p {s,X) for any < X < 1, using a fitting parameter 
a(X). In particular, a(X) acts like a renomarlization 
factor for the number of states of the model and it is 
a(X = 0) = 1 (ferromagnetic case) and a(X) decreasing 
regularly with increasing X . 
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s 


T s 


T v 


7 


6.87 ± 0.04 


6.85 ±0.06 


20 


15.3 ±0.1 


15.3 ±0.1 


50 


31.7 ±0.1 


31.5 ±0.1 



TABLE I. X = 1 and s = 7, 20, 50: Numerical estimates of 
thermodynamic transition temperature T s (s) and percolation 
transition temperature T p (s). For any s they are consistent 
within the numerical error. 
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